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Verfahren und Gerat zur Messung der Biostromverteilung 
HINTERGRTJND DER ERFINDUNG 

Die Erf indung betrifft ein V&rfahren und eine Vorrichtung 
zum Messen der. Biostromverteilung, und insbesondere betrifft 
5 sie ein Verfahren und eine Vorrichtung zum indirekten Messen 
der Biostromverteilung durch Messen des ■ biomagnetischen 
Pelds und durch Abschatzen der Verteilung der Biostrdme, die 
Quellen fur das Magnetfeld sind,. auf Grundlage der gemessen- 
en Daten- 

10 

Die Abschatzung der Biostromverteilung aus den MeSwerten des 
biomagnetischen Felds wird im allgemelnen wie folgt ausge- 
fuhrt. 

15 Sine Beschreibung erfolgt mittels einer neuromagnetischen 
Messung, die eine typ is che Messung fur ein biomagnetisches 
Feld ist. 

Das Koordinatensystem ist so definiert, wie es in Fig. l 
20 dargestellt ist, Es ist angenommen, daB ein a 1b Vektor q n 
reprasentierter Strom an der Position r n in der Zeichhung 
vorhanden 1st. Ein derartiger isolierter Stromvektor kann 
angenommen werden insoweit die Grenzflache des Gehirns auf 
dem Gebiet neuromagnetischer Feldmessung als Kugelf lache be- 
25 trachtet werden kann und dieser Stromvektor wird als Strom- 
dipol bezeiclmet. Das Besugszeichen r TO in Fig. 1 bezeichnet 
die Koordinate des Mefcpunkts . 

Wenn angenommen wird, daB N Strbmdipole im Gehirm existie- 
30 ren, kann der Vektor B m der Magnetf luBdichte am Punkt r m 
mittels des Bio-t-Savart ' s -Gesetzes wie folgt wiedergegeben 
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werden: 



b. = s q ; x(r *- rn) (i) 

n=l 4lT I r,n- r n 13 



Die Magnetfeldkoinponente, die durch ein Magnetometer gemes- 
sen werden kann, ist diejenige Komponente von B m , die recht- 
winklig zur OberflSche des Gehirns steht. Der tatsachlich 
10 bei r m gemessene Wert wird mit D m bezeichnet. 

Wenn.der abgeschatzte Ort und die Stromvektoren jedes Dipols 
dadurch wiedergegeben werden, dafl erhohf das Bezugszeichen 
ahgebracht wird, wird der virtuelle MeBwert D m , der aus 
15 den SchMtzwerten berechnet wird, wie f olgt wiedergebbar : 



Dm = S ( 2 ) . 

n=l 4tt I r ro - f n |3 [ r m l 



20 



Hierbei ist eine Aufwandsfunktion wie folgt definiert: 

?2r f N ; qi, q2f - <3n) 



25 



- S (Dm - D ra )2 v .« . (3) 

m=l • 

Die Gleichung (3) 2eigt das AusinaG der Ubereinstimmung zwi^ 
30 schen einem Schatzwert und einem . tatsachlich gemessenen 
Wert. 

Optiuiale Schatzwerte warden bisher als r lf - - • , r N - g 1# 
q N bestimmt, die die durch Gleichung (3) definierte Auf- 
35 wandsfunktion minimieren. 
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Aus dem Vorstehenden ergibt sich, daS die Anzahl N der bio- 
magnetischen Feldquellen, d.h. der Stromdipole bekannt sein 
sollte, um das herkommliche MeBverfahren ausfuhren- zu kon- 
nen. . Anders gesagt, kann das herkommliche Verfahren nur dann 
ausgefuhrt werden, wenn die Anzahl der Stromdipole bekannt 
ist. 

Wenn die Zahl N, die sich von. der Zahl der tatsachlichen 
Stromdipole unterscheidet, vorgegeben ist und eine Losung 
r l7 . r N ; q lr - q N best immt wird, die die Aufwandsunk- 
tioh minimiert, zeigt die Losung eine Verteilung, die von 
der tatsachlichen Biostromverteilung abweicht- 

ZUS AMMENFAS SUNG DER ERFINDUNG 

ES ist eine Aufgabe der Erfindung, ein verfahren und eine 
Vorrichtung zum Messen der Biostromverteilung zu schaffen, 
die. die Biostromverteilung genau zeigen, und zwar selbst 
dann, wenn die Anzahl der Magnetf eldquellen, d.h. die genaue 
Anzahl der Stromdipole in einem zu messenden Lebewesen nicht 
bekannt ist. Diese Aufgabe wird durcti das Verfahren von 
Anspruch 1 und die Vorrichtung von Anspruch 9 gelost. 

Um die yorstehend angegebene Aufgabe zu losen, werderi gemaB 
der Erfindung mehrere Strome mit einer Anzahl grSBer als 
der jenigen der isolierten Strome, die vermutlich in einem 
Lebewesen vorliegen, angenowmen und es wird ein virtuelles 
Magnetfeld berechnet, das durch den angenommenen Strom am 
0 tatsachlich geraessenen Punkt erzeugt wird. Der angenommene 
Ort und die Stromvektoren des angenommenen Stroms werden ge- 
andert, um die Differenz zwischeh dem virtuellen Magnetfeld 
. und dem aktuellen Magnetfeld zu verringern, um den optimalen 
Schatzwert zu bestimmen. In diesem Fall werden die Absolut- 
5 werte der uberschiissig angenommenen Stromvektoren an null 
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angenahert. 

Wenn die mehreren isolierten Strome/ z. B. Stromdipole mit 
einer Anzahl groBer als der der tatsachlich vorhandenen 
Strome angenommen werden und die Schatzwerte so verahdert 
werden, daB die tfbereinstimmung zwischen der Verteilung des 
aus den angenommenen Stromdipolen berechneten Magnetfelds 
und die Verteilung des tatsachlich gemessenen Magnetfelds 
verbessert 1st, wtirde ein tiberschiissiger Freiheitsgrad e'in- 
gefuhrt. DemgemSB besteht hohe Wahrscheinlichkeit da fur, daB 
eine Losung erhalten wird, die sich von der tatsachlichen 
Stromverteilung unterscheidet . Anders gesagt, ist anzuneh- 
men, daB das am MeBpunkt erzeugte Magnetfeld im wesentlichen 
mit dem durch eineri Stromdipol, der als wahrer stromdipol 
angenommen wird (siehe Gleichung (l) ) erzeugten Magnetfeld 
ttfoereinstimmt, wenn der Strorovektor des Dipols, von dem an- 
genommen wird, daB er an einer Position naher am MeBpunkt 
liegt als es der Position des wahren Stromdipole entspricht, 
kl einer als der wahre Stromdipol ist. 

DemgemaB wird gemaB der Erfindung die Tatsache in. Betracht 
gezogen, daB die angenommenen Stromdipole hinsichtlich der 
Starke die wahren Stromdipole entsprechend der Anzahl' der 
angenommenen Stromdipole gegenuber der Anzahl wahrer Strom- 
vektoren der angenommenen Stromdipole vorliegen, wenri die 
angenommenen Stromdipole nur hinsichtlich der Starke Oder 
der Energie (z. B . des Magnetfelds), zu der sie fiihren, mit 
den wahren Stromdipolen verglicben werden . 

Das heiBt, daB dann, wenn die optimalen Werte der abge- 
schatzten Orte und der Stromvektoren der Strome so bestimmt 
werden, daB die Differenz zwischen dem aus den angenommenen 
Stromen berechneten Magnetfeld und dem tatsachlichen Magnet- 
feld kleiner wird, der uberschiissige Freiheitsgrad der opti- 
malen Werte dadurch verringert wird, daB die Starke der 
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Stromdipole oder die Energie, die sie abgeben, so verringert 
wird, daft die optimalen Werte. der Stromvektoren der angenom- 
. menen Strome, die iiber die Anzahl der wahren Strome hinaus 
vorliegen, auf den Wert null zulaufen. 

5 

Ein uberschUssiger Freiheitsgrad wird dadurch eingeschrankt, 
daS die Gesarotsurnine oder die Gesarotsumme der a-ten Potenz 
der Absolutwerte der stromvektoren der ahgenommenen Stromdi- 
pole oder die Gesaantsunuae der StSrken £er Magnetf elder , die 
10 jeweils alleine durch jeden angenommenen Stroindipol erzeugt: 
werden, oder die Gesamtsumrne der a-ten Potenz der Starken 
derselben yerwendet wird* 



15 



20 



KURZE BESCEREIBUNG DER ZEICHNUNGEN 

Fig. 1 ist eine schematische Ansicht, die ein Ausf utirungs- 
beispiel einer Neuromagnetf eld -Mes sung geiriaB der Erf indung. 
und das. Koordinatensystem veranschaulicht, das zur Abschat- 
zung von Magnetf eldque lien verwendet wird; und 

Fig, 2A und 2B sind FluBdiagramme , die ein Verfahren einer 
simulierten Temperung beirn vorliegenden AusftUirungsbeispiel 
veranscliauliolien . 

25 BESCHREIBUNG DER BEVORZUGTEN AUSFUHRUKGSBEISPIELE 

Es werden nun Ausf Uhrungsbeispiele der Erf indung unter Be- 
zugnahmeauf eine Neuromagnetf eld-Messung beschrieben. Die. 
Neuromagnetf eld-Messung wird dadurch ausgefiihrt, dafc die 

3 0 Komponente des Magnetf elds in einer Richtung rechtwinklig 
zur Gehirrioberf lache an MeBpunkteri urn einen. Kopf herum ge~ 
messen wird, wahrend ein Magnetf luBjnesser 10 von einem Punkt 
12-1 zu einem Punkt 12 -M verstellt wird, wie es in Fig- 1 
dargestellt ist. Die Annaiime - von Raumkoordinaten und Stroro- 

3 5 vektoren von Stromdipolen im Gehirn, die Quellen fur die Er- 
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10 



zeugung von Magnetf eldern sind, kann auf Grundlage der Me-6- 
werte an MeBpunkten um den Kopf . gemafi der Erfindung wie 
folgt geferoffen werden. 

Die Anzahl der MeBpunkte an der Oberf lSche des Kopfs wird 
als M angenommen, und die MeBpunkte werden von l bis m 
durchnumefiert. Die Position des m-ten MeBpunkts ist durch 
den Vektor reprasentiert . Der MeBwert des Magnetfelds an 
diesem Punkt ist als D m reprasentiert. 



Es sei angenommen, da5 N Stromdipole im Him enthalten sind. 
Die Erfindung ist auf einen Fall anwendbar, bei dem die be- 
kannte Anzahl der Dipole nicht .bekannt ist. N wird zeitwei- 
lig alt einer. groBeren Anzahl vorgegeben, als es der Anzahl 
15 von Dipolen entspricht, die verrautlich im Hirn vorliegen. 

Der abgeschatzte ort und die Stromvektoren des so . angenomme- 
nen Stromdipols werden durch r n bzw. q n reprasentiert, wobei 
n = 1, 2, • • * / N gilt. Das Magnetfeld D^ fm = 1, 2, , M) , 

20 das an den MeBpunkten durch die so angenommenen Stromdipole 
erzeugt wird, wird durch eine Berechnung ermitteltV Das 
heiBt, daB als MeSwert eines virtuellen Magnetfelds be- 
zeichnet werden. kann, das durch die angenommenen Stromdipole 
bestimmt wird. Wenn die Oberflache des Kopfs durch eine- Ku- 

25 gelf lache angenahert wird und angenommen wird, daB der Ku- 
gelmittelpunkt der Ursprung der Koordinaten ist, wird D m wie 
folgt mittels einer Gleichung berechnet: 

7 

.30 D m = E — H£_ qti x (r n ,-f n ).^ n , (2J 

n^l 4tt 1 r m -r n |3 I r m \ 

Durch einen Optimierungsvorgang, der im einzelnen beschrie- 
ben wird, werden die abgeschatzten Orts- und stromvektoren 
3 5 der angenommenen Stromdipole so geandert, daB ein Minimal- 
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vert bestimmt wird, den die Aufwandsf unktion einnimmt. 



Beim vorliegenden Ausf uhrungsbeispiel wird die folgende 
Gleichung als Aufwandsf unktion verwendet. 



r 5 N' ^l' q 2' 



2 (D re - D m )2 + Es(q lf q 2 , q H ) (4). 

m=l 



10 



Der erste Term der rechten Forme 1 in der Gleichung (4) be- 
zeichnet die Ubereinstimmung zvischen der Verteilung des 
tatsachlich. gemessenen Magnetfelds und der Verteilung des 
durch die angenommenen Stromdipole erzeugten Magnetfelds. 
15 Wenn die jeweiligen Verteilungen durch Spaltenvektoren G und 
G represent iert warden, kann der erste Term als |G - G| 2 re- 
present iert werden, mit 



20 



G = 



G = 



( I \ 

♦ 



25 Andererseits wird der zweite Term in der Gleichung (4) als 
Gesamtsumme der Magnetfelder angenommen, die jeweils durch 
einen einzelnen angenommenen Stromdipol erzeugt werden. Die 
detaillierte Form, von Es (q x , q 2 q^) i st die folgende: 



30 



Es{q ir q 2r . . , q N ) ~W E I g n 1 

n=l 

wobe i g i 1 1 : 



(5) 



35 



9n = 



(6) 
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und wobei g m 



n wie folgt berechnet wird 



9 m n 



4 * I r n -f n |3 . |r m | 



(7) 



io ... 



Die Komponente g m n in der Richtung' rechtwinklig znr Oberfla- 
che des Gehirns ( Koiripqnente rechtwinklig zur Flache einer 
,15 MeBspule) des Magnetfelds, wie.es am m-ten Mefcpunkt durch 
den n-ten Stromdipol erzeugt wird # wird fur jeden Wert m = 
1, 2, . M durch die Gleichung ( 7 )- bar echnet - Die Starke 
des Spaltenvektors g n mit den jeweiligen Komponenten, d.h. 
.... - der -Absolutwert. |g n | des- durch den -n-ten Stromdipol erzeug- - 

2 0 ten Magnetfelds wird ber echnet,. Genau gesagt, wird die Wur- 

zel der Summe der -Quag.ra.te der Normal)comppriente des durch 
den n-ten Stromdipol- an jedem MeBpunkt erzeugten Magnetfelds 
berechnet. Das Magnetfeld |g n | t wie es individuell durch je- 
den angenommenen Stromdipol erseugt wird, wird auf solche 
25 Weise bestimmt, und die Gesamtsumme hinsichtlich der n ange- 
nommenen Stromdipole wird so bestimmt, da3 sich Es in der 
.Gleichurig (5) ergibt. 

Bs (Ql/ •--/ qjj) f durch die Gleichung (5) represent iert,- 

3 0 dient dazu, den iiberschussigen Freiheitsgrad zu unterdruk- 

ken, wie er dadurch eingefuhrt wurde, daB mehr Dipole ange- 
nommen wurden, als tatsachlich solche vorhanden sind. Obwohl 
es schwierig ist, dies im Fall eines nichtlinearen Umkehr- 
• problems wie bei der Abschatzung der Magnetf eldguelle im Ge- 
35 hirn,.was eine Aufgabe der Erfindung ist, logisch und genau 
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zu verifizieren, kann ein lineares Umkehrproblem matheraa- 
tisch verifiziert werden. Diese Verif izierung ist in "An 
analysis of cost function used in simulated annealing for CT 
image reconstruction", Applied Optics im Druck, von H. Ha- 
5 neishi et al erortert- 

«w" in der Gleichung (5) ist eine Kpnstante zum Gewichten 
des ersten und zweiten Terms im rechten Ausdruck der Glei- 
chung (4) . "w" wird dadurch vorgegeben, daB vorab AE so als 
10 versuchswert berechnet wird, dafi der. Beitrag nicht nur auf , 
dem ersten und zweiten Term in der Gleichung (.4) beruht, 
wenn eine Anderung in der Aufwandsfunktion AE fiir einen Op- 
timierungsvorgang berechnet wird, der nachfolgend beschtie- 
ben wird. 

15 

Obwobl die Gleichung ( 5 ) unmittelbar als Term von Es in der 
Gleichung (4) verwendet werden kann, ist es bevorzugter, ei- 
ne umgeformte Gleichung wie folgt zu verwenden: 

20 N 

Es(q,, q 2 , %)= .« 2 SSnT ($') 

n=l 



25 



wobei der Wert der Potenz a eine positive reelle Zahl ist. 
Die Ro lie von a in der Gleichung (5') ist die folgende: 



Die wahre Stromverteilung wird gemaB der Erfindung dadurch 
abgeschatzt, daB stromdipole angenommen werden, deren Anzahl 

3 0 grdBer als die Anzahl von stromdipoleri ist, die vermutlich 
im Him vorhanden sind. DemgemS/J besteht die Moglichkeit, 
daB gemaB der Abschatzung mehrere Stromdipole im wesentli- 
chen an derselben Position liegen. a in der Gleichung (5') 
. bestimmt die Art der in einem solchen Fall erhaltenen L6- 

35 sung. Es wird. nun beispielhaft ein Fall beschrieben, bei dem 
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zwei Stromdipole u und v im wesentlichen an derselben Posi- 
tion angenommen werden, d.h. ein Fall, bei dem gilt r u == r v . 

(1) Im Fall o < a < 1 wird einer dear Stromvektoren -q u und q v 
5 fiir den u-ten und den v-ten Stromdipol sehr klein. Das 

heiBt, daB |q u | = 0 oder |q v |- =? 0 gilt. Anders gesagt . mehrt 
sioh ein Stromvektor unter mehreren Stramvektoren mit uber- 
lagerten Raumkoordinaten dem Stromvektor eines Stromdipols 
an, der von Natur aus an dieser Position vorliegt, und die 
10 StSrken der Stromvektoren der anderen Stromdipole werden im 
wesentlichen null. 

(2) Im Fall . oe > 1 gilt q u = q v - Dtese Werte entsprechen uri- 
gefahr der Halfte des Werts des stromvektors des Stromdi- 

15 pols, der von Natur aus an der Position r u f r v vorliegt. 

(3) Im Fall a. = 1 itiehrt sich der Wert q u 4* q v dem Wert des 
Stromvektors des Stromdipols an der Position r u = r v an. 
Es kann nicht vorhergesagt werden, welchen Wert ein Verhalt- 

20 nis |q u | zu |q v j annimt. 

Kurz gesagt, ist es bevorzugt, daB <x von 1 abweicht, wobei 
0 < a. < l am bevorzugtesten ist (z. B. a = 1/2) . 

25 Es wird nun im einzelnen ein • Optimierungsvorgang beschrie- 
ben, der SchStzwerte r n , q n (n * 1, N) bestimmt, die 

die durch Gleichung (4) definierte Aufwandsf unktion E (r^, . . 

q N ) , genauer gesagt, die Gleichung (5) oder 
(5 1 ), minimal machen. Die Aufwandsf unktion E ist hinsicht- 

30 lich der Variablen nichtlinear und verftigt iiber einen ortli- 
chen Minimalwert. Demgemafi ist es umnoglich, durch bekannte 
Vorgehensweise fiir nichtlineare Optimierung den wahren Mini- 
malwert (globalen Minimalwert) zu bestimmen. Das vorliegende 
Ausfiihrungsbeispiel verwendet ein in jUngerer Zeit vorge- 

3 5 schlagenes Temperungs-Simulationsverf ahren, das selbst in 
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einera solchen Fall .den globalen Minimalwert bestinunen kann. 

Dieses Verfahren wird nun unter Bezugnahirte auf die Fig- 2A 
und B bescbxieben. 

Die oben angegebenen N Stromdipole werden angenommen und es 
wird ein geeigneter Wert als Anf angswert fiir die jeweilige 
Koordinate r n (n = I, . . / N) vorgegeben, und der Stromvek- . 

tor q n (n = 1, /• N) wird als Variable X n ( n = 1, 2,-.., 

N) vorgegeben, und der Stromvektor q n (n = 1, 2, N) je- 

des Stroicidipols wird als Variable X n (n = N + 1, N + 2/ • • 
2N) vorgegeben. Die Werte der Aufwandsfunktion E f die im 
einzelnen beschrieben wurde, werden unter Verwendung; der 
tatsaciilich. gemessenen und auf gezeichneten Werte D m (m = 1, 
. M) des Magnetfelds berechnet (Schritt 101)- 

Danach wird der Wert von X n um ein sehr kleines Mali geandert 
(dies wird als Versuch bezeichnet) . Unter Bezugnahme auf ei- 
he Anderung AE"der Auf wands funkt ion .Ey wie- sie —dur-eh- den 
Ver such hervorgeru fen wird, wird ermittelt, bb der Versuch 
akzeptierbar ist oder nicht- Jeder Wert X n wird allTnHhlicli 
an den wahren Wert angenahert. 

Als erstes wird ein ausreichend groBer Parameter T vorgege- 
ben/ Dieser Parameter T bestimmt die Wahrscheinlichkeit , daB 
die Xnderung selbst bei AE > 0. akzeptiert wird und er wird 
als Teroperatur bezeichnet. 

AE wird durch eine Versuchsberechnung berechnet und. T wird 
so vorgegeben, daB exp ( -AE / T) = 0,8. bis 0,9 gilt. Die 
Anderungen | Ar | und |Aq| bei einem Versuch werden jeweils 
fiir Variablen r x , - r N und q x , q 2 / <3n vorgegeben 

(Schritt 102) . 

Es ist zweckdienlich, daB der Wert, von | Ar | ungefahr 1 mm 
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' bis 10 mm betragt and der Wert |Aq| ungefahr 1/10 bis 1/1000 
eines vorgegebenen Maximalwerts von jqj betragt. 

Dann werden jeweilige Anfangswerte null fiir Parameter N T , 
5 N Xf N 2 fiir eine wiederholte Einstellung des Versuchs vorge- 
geben (Schritt 103) . 

Eine Variable X n , die Gegenstand des Versuchs ist, wird dann 
ausgewahlt (Schritt 106), um eine Anderung des Vektors AX n 
10 zu bestimmen (Schritt 107). Wenn die Variable .X n eine Koo'r- 
dinate ist, werden die folgenden Gleichungen verwendet. 

(Ax n )x = I if I sin 6 cos ip ^ . 
15 ( Ax n ) y = I A f I sin 9 cos i|i p>. .^gj 

(Ax n ) z = I A r I cos 0 

J 

Wenn die Variable X n ein Stromvektor ist, werden die folgen- 
~den Gieichungetiverwendetv ; — 

20 ~N 

(AXnh = I -A q I sin 0 cos ^ 

(A x n ) Y = I A g I sin 9 cos V (9) 

. ( 

(Ax n )z = Uql cos e 

Die Vektoranderung wird dadurch bestimmt, da3 gleichmaBige 
Zufallszahlen e und ^ in den Bereich o. < e < 2tt bzw. . o < 
* < 7T erzeugt werden . 

30 Die Variable x n wird im Schritt 108 durch X n + Ax n ersetzt, 
Anders gesagt, wird ,ein Versuch vorgenoitiinen . 



Im Schritt 10 9 wird die Differenz AE zvischen der vorange- 
henden und der folgenden Auf wands function E, was zu einer 
35 Anderung AX n fiihrt, wie folgt berechnet: 
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AE = E (x x , . . . , X n + Ax^ , . . . , x 2 h ) 

- E(x lf . x nf . x 2N ) (10) 

5 Wenn AE < 0 gilt, wird die Xnderung im schritt ill ange- 

nommen. Das heiBt, daB das Berechnungsprogramm zum Schritt 
106 zuriickkehrt, um einen Versuch mit der nachsten Variablen 
X n+1 'in Schritten 112 und 105 auszufuhren, wahrend X n + Ax n 
fur X n verwendet wird. Wenn AE > 0 gilt, wird entsprechend 

10 der Wahrscheinlichkeit P(AE) « exp (— AE/T) die von der .Tem- 
peiratur abhangt, bestimmt, ob die Anderung AX n angenommen 
oder zuriickgewiesen wird/ wie es in den Schritten 113, 114 
und lis veranschaulicht ist. Wenn bestimmt wird, daft , die . Xn- 
. derung angenommen wird, kehrt das- Berechnungsprogramm uber 

15 die Schritte 115 und 105 zum Scliritt 106 zuriick, wahrend 
X n + AX n als neuer Wert'X n verwendet wird. Wenn bestimmt 
wird, daB die Anderung zuriickgewiesen wird, werden die Auf- 
wands funkt ion E und die Ur sprung swerte im Schritt 117 wie- 
derhergestellt und das Programm kehrt iiber die Schritte 118 

20 und 105 zum Schritt 106 zuriick. 

Eine derartige Fo.lge von Schritten wird fur jede Variable 
ausgefiihrt, d.h- fur n = 1 bis 2N. Dieser Vorgang wird ins- 
gesamt N T max Male wiederholt. Es ist zweckdienlich, daB z. 
25 B. ein Wert von ungefahr 100 bis 400 flir N T TOax verwendet 
wird. 

im Schritt 112 des FluBdiagramms von Fig. 2 A bezeichnet 
die Anzahl von Anderungen, die fur diejenige Richtung ange- 

30 nommen werden, bei der sich die- Auf wands funkt ion verringert, 
und N 2 im Schritt 115 bezeichnet die Anzahl von Anderungen, 
die fur diejenige Richtung angenonunen . wurden, fur die die 
Aufwands funkt ion zunimmt. N3 im Schritt 105 bezeichnet die 
Anzahl zuruckgewiesener Anderungen. Wenn sich im Schritt 105 

35 herausstellt, daB N T max Versuche abgeschlossen sind, geht 
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das Programm 2u einem Schritt weiter r in dem die Ermittlung 
ausgefiihrt wird, ob die die Annahmewahrscheinlichkeit fur 
eine Anderung bestimmende Tempera tur T verringert wird« . Das 
heifct, daB auf Grundlage des Werts \N ± - N 2 1 bestimmt wird, 
5 ob die Anzahl N L von Anstiegen der Auf wands function bei den 
angenoramenen Versuchen naherungsweise der Anzahl von Verriri-. 
gerungen der Auf wandsf unktion entspricht. Wenn dieser Wert 
nicht kleiner als e (z. B. ungefahr 0,02) ist, kehrt das 
Programm zum Schritt 103 zuruck, ohne die Temper atur zu an- 

10 dem, und Versuche werden N-r 10 ^ Male wiederholt. Wenn der 
Wert |Nx - N 2 | / N kleiner als e ist, geht.das Programm zu 
einem Schritt 120 weiter. Durch Verringern der Temperatur ' T 
wirdi <iie Wahrscheirilichkeit-, daB eine Anderung zum Erhalten 
von AE > 0 angenommen wird, verringert, und danach.geht das 

15. Programm zum Schritt 103 weiter und die Versuche werden auf 
ahnliche Weise N T TOa * Male wiederholt. 

Beim vorliegenden Ausfuhrungsbeispiel wird zum Absehken der 
Temperatur im Schritt 12 0 eine Vorgehensweise verwendet, die 

20 versucht, im k-ten Temperaturschritt T^ — fT^^ zu erhalten. 
Dabei wird fur £ ein Wert von ungefahr 0,9 bis 0,95 verwen- 
det. In einer Literaturetelle (H. Szu et al r "Fast simulated 
Annealing", Physics Letters A, Vol. 122, S. 157, 1987) ist 
auch vorgeschlagen, die Temperatur urn T^ = T 0 / (1 + k) oder 

25 T k = T 0 / log (e + k) zu verringern. 

Da N T max Jinderungs versuche fixr jede Variable, wie sie mit 
einigen Zyklen wiederholt werden, allmahlich die Temperatur 
T verringern, wie oben angegeben, nahert sich jede 'variable 

3 0 einem optimalen Wert an und die angenommene Jiriderung nimmt 
allmahlich ab. Wenn im Schritt 118 + N 2 = 0 bestimmt 

wird, d.h. wenn keine angenommene Anderung bei den versuch- 
ten Anderungen besteht, wird eine Variable erzielt, die die 
Auf wandsf unktion minimiert, d.h. optimale Werte fur den Ort 

35 und die Stromvektoren eines Stromdipols. Daher wird der si- 
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roulierte T^mperungsvorgang beendet. 

Obwohl das vorliegende Ausf tthrungsbeispiel unter Bezugnahme 
auf einen Fall beschrieben wurde, bei dem simulierte Tempe- 
rung als Optimierungsverf ahren verwendet ist, 1st die Er fin- 
dung nicht auf diesen Fall beschrankt. Andere Algorithmen, 
die einen wahren Minimalwert einer Punktion auffinden kon- 
nen, die. einen ort lichen Minimalwert entfcalt, konnen verwen- 
det werden, wie der von J. K* Holland vorgeschlagene gene- 
tische Algoritixmus . 

Wenn ' Bere.iche, in denen Dipole existieren, vorab bekannt 
sind, kann diese A-Pr iori-Xnf orjnation wie folgt in den Opti- 
xnierungsvorgang eingebaut werden. 

in dieserc Fall wird jeder Schatzwert/ der E minimiert, durch 
die folgende Gleichung (11) anstatt durch die Gleichung (4) 
bestimmt. 

E(r ir r 2 r .-r %L' '2' *' ^ 

m=l 

+ E L ( £ 1' f N } (L1) 

El 1st wie folgt definiert: 

1) Wenn alle Werte r x , . r N in vorab eingestellten Berei- 
chen existieren, gilt E L (f lf . - r N ) *= 0. 

.2) Wenn einer der Werte x lt r N aus. den vorab einge- 

stellten Bereichen herauslauft, gilt E L (r lf . r N >- = A,, 

wobei es sich ura einen Wert handelt, der fur die aktuelle 
Temperatur T 1 ausreichend hoch ist, d.h. einen Wert, bei dera 
gilt exp (-A / T) =0. Wenn E L Auf. diese Weise definiert 
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wird, ist die Wahrscheinlichkeit , da£ eine Anderung angenoin- 
roe^n^/ird, bei der jedes r n aus dem vorab ... eingestellten Be- 
reich herauslauft, ungefahr null. 1m vorab eingestellten Be- 
reich kann ein Schatzwert ermittelt werden, der dexr folgen- 
5 den Wert miniiaiert : 

B(£ ir . . , f N ; r . . . t q M ) 

10 = S (Q» - D m )2 + Es(q x/ - • ♦ ,.q H ) 

m-1 

Obwohl ein Fall beschrieben wurde, bei dem der durch die 
Gleichung (5) oder (5 r .) ausgedriickte Wert, d-h. das durch 
den angenommenen Stromdipol erzeugte Magnetfeld als * zweiter 

15 Term als Gleichung (4) verwendet wird, urn den Preiheitsgrad 
zu beschrSnken, kann anstelle der Gleichung (5} oder (5 % ) 
eine Gleichung (5 1 » ) entsprechend der Gesairitsumme der Abso- 
lutwerte der Stromvektoren der angenomitienen Stromdipole oder 
eine Gleichung (5M.«) der u-ten Potenz der Gesamtsumme ■ ver- 

2 0 wendet werden: 

K 

EsCq^ q 2 , , q N ) = w S I q n I (5") 

n=l 

25 n 

Es{§ ir q 2 , . q H ) = w 2 I q n (a (5 '") 

n=l 

In den. Gleichungen (5' 1 ) und (5 ,f ') sind die Bedeutungen von 
w und a dieselben wie in den Gleichungen (5) und (5 1 ) . 

30 

Da Terme von Positionsvektoren in den Gleichungen (5) und 
(5 1 ) enthalten sind, ist der Preiheitsgrad starker be- 
schrankt als in den Gleichungen (5 1 9 ) und <5 1 ' ' ) . 
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PATENT ANS PROCHK 

1. Verfahren zur Messung einer Biostrom-Verteilung. durch 
Messen des von der Aktivitat eines Lebewesens erzeugten bio- 
magnetischen Feldes an einer. Vielzahl von Mefipunkten (12- 
1,...,12-M) und Bestiinmen isolierter Strome (c^) in dem Lebe^ 
wesen aus den gemessenen Werten des Hagnetfelds, mit folgen- 
den Schritten: 

Messen des tatsachlichen Magnetfelds an den genannten 
Meftpunkten, 

.Annehmen einer Vielzahl von Stromen in einer gr5£eren 
Zahl als deir erwartungsgemaB in dem Lebewesen vorhahdenen 
.Zahl isolierter Strome und Vorgeben geschatzter Orte und 
Stromvektoren der angenommenen Str5iae, 

Berechnen eines virtuelleii Magnetf elds an jedem Meftpunkt 
aus den geschatzten Or ten und Stromvektoren, 

Berechnen einer Kostenfunktion, die die Summe eines 

erst en und eines zweiten Terms bel nha ltet, von denen der 

erste Term den Unterschied zwischen den virtuellen und tat- 
sachlichem Magnetf eldern darstellt und der zweite Term aus 
den geschatzten Stroiavektoren oder den von diesen jeweils er- 
zeugten Magnetf^ldem berechnet wird und den durch die Annah- 
iae von mehr Stromen, als . in dem Lebewesen erwartungsgem^B an 
isolierten Stromen vorhanden sind, eingefuhrten Freibeitsgrad 
unterdriickt, und 

Miniiaiereii der Kostenfunktion durch fortlauf endes Andern 
jeder der geschatzten Orte und, Stromvektoren der angenommenen 
Strome zur Annaherung an die Orte und Stroiavektoren der in 
dem Lebewesen vorhandenen isolierten StrOme. 

2. Verfahren nach Anspruch 1, wobei der zweite Term die Ge- 
samtsumme der Absolutwerte der geschatzten Stroiavektoren dar- 
stellt. 
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3. Verfahren nach Anspruch l f wobei der zweite Term die Ge- 
. samtsumme der a-ten Potenz der. Absolutwerte der geschatzten 

Stromvektoren darstellt und oc.eine positive Realzahl 1st. 

4. Verfahren nach Anspruch 1, wobei der zweite Term die Ge- 
samtsumme der ■ Absolutwerte der an jedem MeBpunkt von jedem 
der angenommenen Stromvektoren erzeugten Magnetf elder dar- 
stellt. 

5.. Verfahren nach Anspruch 1, wobei der zweite Term die Ge~ 
samtsuinme der Absolutwerte der senkrecht zur OberflSche des 
Lebewesens gerichteten Magnetf eldkomponenten der von jedem 
. der angenommen Stromvektoren an jedem Meflpunkt erzeugten 
Magnetf elder darstellt, 

6„ Verfahren nach Anspruch 1, wobei der zweite Term die 
Gesamtsumme der a-ten Potehz der Absolutwerte der senkrecht 
zur Oberf lS.che des Lebewesens gerichteten 

Magnetteldkomponenten der von jedem der angenommenen 
Stromvektoren an jedem MeBpunkt erzeugten Magnetf elder . 
darstellt und a eine positive Realzahl 1st. 

7. .Verfahren nach Anspruch 3 oder 6 f mit a * 1, 
vorzugsweise 0 < a < 1. . 

8, Verfahren nach einem der ; - Ansprtiche 1 bis l r wobei der 
Schritt zur Minimierung der Kostenfunktion einen ersten 
Unterschritt zur Annahme einer Anderung der geschatzten Orte 
und Stromvektoren dann, wenn sich die Kostenfunktion 
verringert, und eineh zweiten Unterschritt zur Annahme einer 
Anderung der geschatzten Orte und Stromvektoren mit eiher 
Wahrscheinlichkeit von exp(-AE/T), wobei AE eine Anderung 
der Kostenfunktion und T einen vorbestimmten Parameter, der 
die Wahrscheinlichkeit einer Annahme der Anderung. festlegt, 
auch wenn sich die Kostenfunktion erhoht, aufweist. . 
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9. vorrichtung zur Messung einer Biostromverteilung. durch 
Messen des von der Aktivitat eines Lebewesens erzeugten 
magnetischen Felds an einer Vielzahl von Mefipunkten (12- 
l...,12-M) und Bestimmen der Verteilung isolierter Str5me 
(qn) in dem Lebewesen aus den gemessenen Werfcen des 
Magnetfelds, aufweisend: 

eine Einrichtung zur Messung des tatsachlichen 
Magnetfelds an den genannten MeBpunkten, 

.eine Einrichtung zur Annahme einer Vielzahl von StrOmen 
in einer grGGeren Zahl als der Anzahl erwartungsgemaft in dem 
Lebewesen vorhandener isolierter Strome und Vorgeben 
geschatzter Orte und Stromvektoren der angenommen Strome, . 

eine Einrichtung zur Berechnung eines virtue-lien 
Magnetfelds an jedem der MeJSpunkte aus den geschatzten Orten 

und Stromvektoren, 

eine Einrichtung zur Berechnung einer Kostenf unktion die 
eine Summe aus einem ersten und einem zweiten Term 
beinhaltet , von denon der erste Term die Differenz zwischen - 
dem virtuellen und dem tatsachlichen Magnetfeld darstellt und 
der zweite Term aus den geschatzten Stromvektoren Oder den 
von diesen jeweils erzeugten Magnetf eldern berechnet wird und 
der Unterdruckung des dadurch eingefiihrten Freiheitsgrads 
dient, dafi mehr Strome angenommen wurden als erwartungsgemafi 
an isolierten Stromen in dem Lebewesen bestehen, und 

eine Einrichtung zur Minimierung der Kostenf unktion 
durch fortlaufendes Andern jeder der geschatzten Orte und 
Stromvektoren der angenommen Strome zur "Annaherung an die 
Orte und stromvektoren der isolierten Strome in dem 
Lebewesen. 
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